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Abstract 

We study the partition function of a two-dimensional Coulomb gas on a 
circle, in the presence of external pointlike charges, in a double scaling limit 
where both the external charges and the number of gas particles are large. Our 
original motivation comes from studying amplitudes for multi-string emission 
from a decaying D-brane in the high-energy limit. We analyze the scaling 
limit of the partition function and calculate explicit results. We also consider 
applications to random matrix theory. The partition functions can be related 
to random scattering, or to weights of lattice paths in certain growth models. 
In particular, we consider the discrete polynuclear growth model and use our 
results to compute the cumulative probability density for the height of long 
level- 1 paths. We also obtain an estimate for an almost certain maximum 
height. 
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1 Introduction 

Unstable D-branes have many decay channels. The simplest, and a relatively well 
understood one, is closed string emission, one at a time at tree level |l}j4]- The multi- 
closed string or multi- open string production channels lead to more complicated 
calculations, and are less well understood. In a sequence of papers (5j [6] we have 
developed estimates for amplitudes in the multi- open string channels, and most 
recently in the closed string pair production channel [?]. Our calculational strategy 
was based on mapping the problem to calculating the partition function of a log-gas 
of unit charges in a unit circle (the Dyson gas [8], arising from the background of 
the condensing tachyons on the brane) with additional test charges (arising from the 
vertex operators for string emission) [9-11 . We then used electrostatics to find the 
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leading saddle-point contribution. This is sufficient for the high-energy limit of the 
emission amplitude [7]. 

Logarithmic Coulomb interaction potential arises as a solution of the Poisson equa- 
tion in two dimensions. After the original paper [8 , thermodynamics of more general 
Coulomb systems have been studied, e.g., in 12 . Some applications of Coulomb 
gas thermodynamics are discussed in 13 . Furthermore, Dyson gas with additional 



charges has applications in various contexts in random matrix theory (RMT), such 
as estimates for scattering amplitudes in a quantum mechanical random scattering 
problem, and estimates for the weights of long lattice paths in growth models. It 
is therefore worthwhile to study if the electrostatics approach can be used to derive 
explicit results. 

In this paper we analyze the electrostatics approach in more detail. We will 
also discuss applications to random scattering and to growth models. As a first 
example we consider the boundary one-point partition function, with one test charge 
on the boundary. This configuration is solved exactly in the literature, and serves 
as a test for the electrostatic method. We point out that the one-point function 



is related to lattice paths in the "lock step model of vicious walkers" 14-16 . We 
also derive new results for higher point boundary functions. Their application to 
open string emission from decaying D-branes will be discussed elsewhere. We then 
consider bulk n-point functions with n external charges inside the unit circle. In string 
theory, they arise from (multi) closed string emission from decaying D-branes. A more 
widely known physics application is multi-channel scattering across a non-ideal lead 
(with a tunneling barrier) coupled to a chaotic cavity. We show that the bulk one- 
point partition function, which was analyzed in detail in our earlier work [7], can be 
related to the Poisson kernel associated with the eigenvalue probability distribution 
function for the scattering matrix. Another RMT interpretation, applicable for n- 
point functions for all n > 1, is given by weights of a class of long lattice paths 
in the discrete polynuclear growth (PNG) model (we follow the presentation in |l6j 
1~7]). Using the results from our electrostatic analysis in |7j, we give an estimate for 
the cumulative probability distribution function of long paths from large weighted 
matrices. This limit also corresponds to the high-energy limit of closed string pair 
production from a decaying D-brane [7|. 

As the above examples hopefully show, there is an interesting connection between 
D-brane decay in string theory and problems in random matrix theory. We antic- 
ipate that there are many other examples to be found. The connections help in 
finding useful strategies for problem-solving in string theory. They should also serve 
as additional motivation for work in random matrix theory. For example, in the 
string theory context one is naturally lead to consider log-gases in grand canonical 
ensembles, whereas our impression is that most of mathematical work is focused on 
canonical ensembles. A very recent example of potentially related work on grand 



canonical ensembles is 18 



The sectioning of the paper is as follows. In Section [2] we introduce our approach 
and show that the leading behavior of the partition functions can be calculated by 
using continuum electrostatics. We also analyze the associated electrostatic potential 
problem. In Section |3j as a first example we calculate in detail the leading behavior 
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of the partition function with one external charge on the unit circle, and compare 
with the available exact result. We also point out the lattice path interpretation. In 
Section [4] we move to consider boundary n-point functions, analyzing the two-point 
function in detail. In Section [5] we move to consider external charges in the bulk, 
starting with the one-point function. Depending on the chosen parameters, this case 
is either very simple, with application to random scattering, or more complicated but 
where we can use the electrostatic approximation in the double scaling limit. Finally 
in Section [6] we consider bulk n-point functions and their relation to the discrete 
polynuclear growth model, and calculate the cumulative probability density function 
in a special case. We also discuss a connection to lattice paths. Appendix |A| contains 



calculational details. In 19 we apply the results of Section 4 to D-brane decay. 



2 Large N expansion of the partition function 

We consider the partition function^ of a two-dimensional Coulomb gas on a circle 
(Dyson gas) with n external charges £ a at fixed positions w a on the complex plane 5][(r 



N 

I I- 



' J i=l 



where the Hamiltonian H reads 

N n 

H = - log|e^-e^|-^^^log|e*- Wa |- ^ Ublog\w a -w b \ . (2) 

l<i<j<N *=1 a=1 l<a<b<n 

It describes an ensemble of iV unit charges on the unit circle at e*' 1 , . . . , e itN , which 
interact via the Coulomb potential and are subject to the electric field created by the 
fixed test charges £ a . The charges are located at w a anywhere on the complex plane. 
The charges with \w a \ = 1 {\w a \ < 1) will be called boundary (bulk) charges. For the 
boundary charges we also denote w a = e lTa . 

A useful strategy is to study the system using classical electrostatics, first in the 
limit of a large number A" of unit charges. In this limit the unit charges can be 
treated as a continuous charge distribution. Its density is modified by the external 
charges £ a , and we are interested in finding the electrostatic equilibrium configuration. 
In particular, nontrivial electrostatic configurations are obtained when the external 
charges are of the same order as the sum of the unit charges, A" ~ £ a . We shall in 
fact study the behavior of the system in a double scaling limit limit N, £ a — > oo with 
all ratios N/£ a fixedQ 

The partition function in this limit (at (5 = 2) can be used to analyze the emis- 



sion of highly energetic strings from a decaying D-brane. This is reminiscent of 21 



lr The full partition function factorizes to Z = Zki n Z pot where the contribution from kinetic 
energies Zkm involves only simple Gaussian integrals. Hence we focus on the non-trivial potential 
energy contribution only. 

2 The behavior of Z n in the limit N — > oo while holding £ a fixed can be calculated using the 
asymptotics of Toeplitz determinants [9 20 . 
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where Gross and Mende used electrostatic equilibrium conditions to identify the dom- 
inant saddle-point contributions to string scattering amplitudes at arbitrary order to 
determine their high-energy behavior. 

The system is analogous to the circular unitary ensemble (CUE) of random matrix 
theory (an early observation of this is [22]), where the connection to the Coulomb 
gas is a standard tool. Our method indeed resembles the electrostatic derivation of 
Wigner's semi-circle law [23] (see also [16]). The partition function ([T| is related to 
an expectation value of a periodic function in the circular unitary ensemble (CUE), 

n N 

iti \0€a 



Z n ((3; N) = U K - w b \^ ( n n I 1 



a<b \a=l i=l 




where /3 = 1,2, 4 correspond to ensembles of symmetric, unitary, and self-dual quater- 
nion unitary matrices U(N)/0(N), U(N), and U(2N)/Sp(2N), other values define 
the general /3-ensemble|^] 

Let us expand the partition function at large N, with N/^ a fixed. We anticipate 
the following series expansion: 

\ogZ n {(3; {&}; N) = C N 2 + C,N l + C 2 N° + • ■ ■ , (4) 

where the coefficients Ck generally depend on the ratios N/^ a and may contain log- 
arithms of N. Using a saddle-point-like approach we divide the exponent into the 
leading contribution and fluctuations, 

logZ n (/?;{U;A0 = -PH + A f 

= -P£ + A d + A f , (5) 

where Ho is the minimum value of the Hamiltonian ^ and Af arises from "fluctu- 
ations" around the minimum value of the Hamiltonian. We further divided Hq into 
£, the (total) electrostatic energy of the minimum configuration approximated with 
continuous charge densities, and the corrections from the difference of the energies 
of the discrete and continuous systems. These corrections are analyzed in detail in 
Appendix |A] and they are shown to contribute only at next-to-leading order [O (N)]. 
Hence the leading term CqN 2 is determined by the first term — (38 in Indeed, the 
energy £ contributes only at O (N 2 ), so 

C N 2 = -f3£ , (6) 

and the next-to- leading coefficient C\ arises from the corrections A (see Appendix [A|. 

2.1 Calculation of the leading term 

We now discuss how the leading term CqN 2 = —f3£ is calculated. The first ingredient 
is the (complex) potential V(w). The physical potential U is the real part, 

U{w) = ReV(w) . (7) 



3 Random matrix interpretations of the generic circular unitary /3-ensembles have been constructed 
in [24l. 
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The potential is sourced by the continuous charge density pit) (which we use to 
approximate the N unit charges in the N — > oo limit, with the normalization iV = 
§ dtp{t)) and the external charges £ a , 



V(w) = - J dtpit) log (w - e lt ) - ^f a logl 



w — w„ 



We fixed the zero level of the potential by setting the possible constant term to zero. 
We then need to find V(w) and p(t) in the equilibrium configuration. The continuous 
density consists only of positive elementary charges. 

A characteristic feature of the ensuing electrostatic configurations is the presence 
of gaps in the charge distribution p(t). Since the distribution is made of equal (posi- 
tive) unit charges, p(t) cannot become negative. Hence, the charge distribution may 
vanish in the vicinity of an external charge if it is close enough to the unit circle. A 
detailed discussion of this phenomenon in the case of bulk two-point amplitude can be 
found in [7]. Thus, we are interested in configurations where pit) is strictly positive 
within a domain C of the unit circle and vanishes elsewhere; the gaps are created by 
the repulsive force from the positive external charges £ Q . 

The problem of finding the equilibrium configuration becomes a standard poten- 
tial problem of classical electrostatics in the presence of external charges £ a and a 
conductor which fills the domain C . In the equilibrium configuration the physical 
potential U (e lt ) must be constant within C, 

U(e a ) = U , e lt eC. (9) 

If the potential would vary within C, the charges on the conductor would be subject 
to forces, and move until the potential adjusts to a constant value. 

The electrostatic energy is defined as the energy of the system in the equilibrium 
configuration, and may be expressed in terms of po(t) as 

£ = ~Jdt 1 dt2Po{ti)po(t2)los\e itl -e it '\ 

n „ 

~y^Aa dtp (t) log \e %t - W a \ - ^ log \w a ~ W b \ 

a=l J l<a<b<n 

= £ c + £^ + £ ( . (10) 

We again set a possible constant contribution to zero so that the energy matches 
with the continuum limit of the Hamiltonian ^ . By using the fact that the physical 
potential U{e %t ) is constant within the domain C of nonzero po(t), we can express the 
result as 

N 1 

£ = jU + -£ c ^ + £s, (11) 

where the first term includes one half of the "interaction energy" £ c -£- Hence, to 
calculate £ and the leading behavior of Z n for large iV ~ £ , it is sufficient to find Uq 



and £ c _£. We present a general method for this in Section 2.2.1 In the rest of this 
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paper we will drop the subscript from the charge density and the potential, so p(t) 
and U will refer to the quantities at equilibrium. 

For a system with several external charges many distinct gaps are possible. Their 
locations are fixed by the requirement that p(t) must vanish continuously at their 
endpointsj^] We will next present the potential problem and find explicit solutions. 
The reader who is interested in applications of the solutions may wish to jump to 
Section [3] upon first reading. 



2.2 The potential problem and its solution 

We are interested in solving potential problems where the conductors are pieces of 
the arc of the unit circle. We find it useful to map the disk onto the upper half plane 
(UHP) by the conformal map (see Fig. [TJ 



w i— > z = it = q(w) 

1+w v ' 



:i2) 



with the understanding that we work in the compactified complex plane with the 
oo-point included. Here ( is a real parameter. 



Li 



z = i cot( t 12) - 




Figure 1: The conformal map that takes the conducting plane to the real axis. A 
possible mapping of a pointlike charge Q is also illustrated. 



2.2.1 Electric field construction 

Let us first review a powerful approach which allows us to find compact expressions for 
electric fields that correspond to configurations with arbitrary numbers of produced 
strings. The general potential problem has as conductors n separate arcs of the unit 
circle. It is again useful to map the unit circle on the real axis. The analytic properties 
of the complex electric field allow us to write down the solution in a special form [16] . 

4 If the number of conductors h > 2, the domain of nonzero p(t) is disconnected, and the potential 
differences between the connected sub domains are also free parameters. In this case, energy is 
minimized when these differences are all zero. 
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Let us assume for a moment that we have two conductors, n — 2, which lie at [a, b] 
and at [c, d]. We can then define the following functions 

g(z) = yj z — ay/ ' z — by/ z — cyj ' z — d 

h(z) = y/ z — ay/b — zy/c — zy/d — z , (13) 

where the principal branch of the square root function is used. Then g(z) is ana- 
lytic everywhere except on the conductors where it has branch cuts, whereas h(z) is 
analytic over the conductors, but has branch cuts in the gap regions. 

Generalization to n conductors is found by adding more terms to g(z) and h(z). 
We choose the branches such that 

h(z) = Tig{z) , (14) 

where the minus sign holds in the upper half plane and the plus sign holds in the 
lower half plane. It is then possible to show that, when exposed to an appropriate 
external (conjugate) electric field E cxt (z), a solution for the (conjugate) electric field 
due to the conductors reads 

rp f \ 1 f , [ Rggextg) 

E c (z) = -g(z) / dt , (15) 

K Jcond {Z-t)h(t) 

where the integration is over the conductors. Notice that the solution is bounded 
at the conductor endpoints. It is well defined everywhere except on the conductors, 
where it has discontinuities. The charge density is proportional to the discontinuity 
and becomes 

i \ 1 w x f A 1 ReE cxt {z) - ReE ext {t) 

p(z) = ^ h{z) L dt m — — — • (16) 

The 2n endpoints of the conductors are not arbitrary parameters but are subject 
to constraints. In particular, these parameters (the constants a, b, c, and d in the 
h = 2 case) are partially fixed by the asymptotics of the electric field 

E c (z) ~ j (17) 



for z — > oo where N is the total electric charge on the conductors. This gives 

f , t k ReE ext (t) 
= dt V ; , fc = 0,...,n-l 

icond h(t) 
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f , r ReE cxt (t) . . 

icond h-V*) 



There are n — 1 parameters left free, which are identified as the potential differences 
between the conducting planes. The minimum energy is reached when the potential 
differences are equal to zero. Therefore, we require 



Re / dzE(z) = 0, k = 2,...,h, (19) 
Jc k 
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where 

E(z) = E c (z) + E cxt (z) (20) 

is the total electric field and the curve Ck connects the first conductor with the fcth 
one without intersecting the other conductors, for example. 

For the case of interest to us, the external field is due to n point-like charges, 



n 



a=l 



(21) 



For this configuration the integrals above can be calculated exactly by using contour 
integration methods. We find 



E c (z) 
E(z) 



9Kz) 
2 

2 



a 



+ 



g(z a )(z-z a ) g(z*)(z - z* a ) 

1 1 

+ 



Kz) 
2tt 



_g(z a )(z-z a ) g(z*)(z - z* a ) 
1 1 



T<% - Za )( z - zt ) 
^ U (z-z a )(z-z*) 



+ 



g(z a )(z-z a ) g(z*)(z - z* a ) 



(22) 



where z a = x a + zy . The asymptotic conditions imply 



o = 5> 

a 



fe = 0, . . . , n — 1 



(23) 



The potential and the total energy of the system are of special interest to us. We 
can integrate the complex electric field in (|22|) to give 



V c (z) 
V(z) 



[Z — Z n 



a 

a 

\ a / 



+ 



_g( z a)(t-z a ) g(z*){t-z*)_ 



(24) 



Notice that this solution has the asymptotics required by our definition (|8]) (or more 
precisely by its natural generalization to the half plane) 



V(z) = - \N + hgz + O 



(25) 



so that the constant term of the expansion in 1/z is zero. 

Let us then write down explicitly the solutions for the disk, %. e., apply the transfor- 
mation q(w) of ( 12 ) (with, say (3 = 1). To do so, we should start with the configuration 
where charges £1, . . . ,£„ are located at arbitrary points z\ = q(w\), . . . ,z n = q(w n ), 
and an additional charge £ n+ i = — N — J2 a £a * s located at z = q(oo ). Let us denote 
the corresponding electric field and potential, as given by eqs. (22) and (24), by E 
and V, respectively. 

The electric field on the disk may be then written as 



E(w) 
EJw) 



E{q{w)) 



N 



q(w) -q(oo) 



(26) 



and integration gives the potential on the disk 

V(w) 
V c (w) 



V(q(w)) — lim 



dzE(z)+ (iV + VeJlogi? 

J m V a J 



V c (q(w)) + N\og [q(w) - q(oo)\ 

R 



— lim 



dzE(z) - N log R 



q(oo) 



-N\ogj t q(l/t) 



(27) 



t=Q 



Here the integrals are independent of w. They were added to achieve the correct 
asymptotics, 



V(w) 



log w + O 



V c (w) = -Nlogw + O 



(28) 



Finally, the potential on the conductors is found as 

U = Re V(w ) , (29) 
where wo is any point on the conductors, and the "interaction" energy is given by 

4-£ = ^£aReK(0 • (30) 



The total energy then reads (see (11) above) 

£ 



(31) 



2.2.2 Mapping approach 



An alternative and complementary approach is found by mapping conformally the 
complement of the conductors to such a region where the potential solution is known. 

In the case of one conductor, say at argu> = t c . . . 2n — t c , we may map the 
conductor to pieces of the real axis ] — oo, —a] U [a, oof, by choosing ( = acot(t c /2) 
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in (12) (see Fig. [I]). For our purposes it is sufficient to solve for the Green function 
for this configuration, i.e., the solution for one external charge at an arbitrary point. 

At this point we could use the formulas of the previous subsection to directly 
construct the electric field and the potential. However, we may also apply a further 
transformation 



z i — y % 



a + z 



a 



z 



(32) 



that maps the complex plane to the upper half plane and the conductor as the real 
axis. The Green function is then found by inserting a mirror charge on the lower half 
plane and mapping back. The obtained solution reads (for an external charge Q at 

* = z o) 



V(z) = -Q log 



a+z 
a—z 



a+ZQ 



a+z 



+ 



a+zt 



(33) 



where the physical potential is the real part. 

The solution of the potential problem of a special interest to us (with the conductor 



at real axis) is the sum of (33 ) with zq = 0, Q = £ and zq = Q = —N — £, where 



— i£ = — mcot (t c /2) is the image of w = oo in the conformal map (12): 

+ (iV + £)log 



V(z) 



-Clog 



a+z j 

a—z 



a+z 



+ 1 



a+z 


a-K 


a—z * 




a+z j_ 


/a+i( 


~z + 1 


J a-i( 



+ v . 



(34) 



In (34) the Vq is a constant. The charge density is found as the discontinuity of the 



electric field —V'(z)* on the real axis: 

M - 



tc\/z 2 — a* 
£(Wz 2 - a 2 



^2 



z 2 + C 



where 



7ra\z\ (z 2 + ( 2 ) 

C _ y/N(N + 2£) 
a £ 



z\ > a , 



(35) 



(36) 



guarantees that the charge density vanishes at z — ±a. By using the map (12), the 
complete solution with the conductor on the unit circle can be written as 



V{w) 



P(t) 



-Clog 



— w+e ztc 



- 1 



we" c —l 
— w+e lt c 



+ (iv+e)io g 



we it c _l jL-t+r 

»z — e 2 



-w+e" c 



we itc — l 

—w+e %t c 



-Re \ /cot 2 — 
2tt V 2 



cot z 



+ V (37) 



where sin(t c /2) = £/(N + £), and t = aigw. Notice that the charge density vanishes 
for t = —t c . . . t c as the square root is purely imaginary. We use the standard principal 
branch (Re(-y/z) > if | argz| < it) for the square root. 
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2.2.3 Two-gap case 

Let us then go on to the much more complicated case of two conductors. With- 
out any loss of generality, we may assume that the gaps lie at ]— 1/a, l/a[ and at 
]— oo, — a[ U ]a, oof, where a > 1. Hence we need to solve the potential problem for 
two conducting planes at [—a, —1/a] and at [1/a, a]. Boundary values of the potential 
could be taken to be different but we assume that they are equal as this will minimize 
the energy. 

It is useful to further map the problem to a rectangular region. The mapping can 
be constructed by using the Schwarz-Christoffel formula: 

z — y u = p(z) oc / — = oc sa~ 1 (az\m) , (39) 

Jo ^{z> + a)(z' + l/a)(z' - l/a)(z' - a) 

where sn -1 is the inverse of the Jacobi elliptic sine function sn and m = a~ 4 is the 
(elliptic) parameter. We fix the normalization to 



p(z) = sn (az\m) . (40) 



Then we find 



p(0) = 

p(±l/a) = ±K(m)=±K 

p(±a + ie) = ±K(m)+iK(l-m) = ±K + iK' 

p(ioc) = iK(l - m) , (41) 

where K is the complete elliptic integral of the first kind. Hence the upper half plane 
is mapped on the rectangle }—K,K[ x ]0,i^'[ such that the conductors form the 
vertical edges (see Fig. [2]). Since the gap at ]— 1/a, l/a[ is mapped to the real axis, we 
may use the Schwarz reflection principle to argue that the lower half plane is mapped 
to ]—K, K[ x ]—K', 0[ with again conductors on the vertical edges (if the mapping is 
continued analytically through this gap). Interestingly, if Schwarz reflection principle 
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is used to analytically continue the mapping through the other gap, we see that the 
lower half plane is also mapped to ]— K, K[ x ]K',2K'[. Repeating the analytic 
continuation through the gaps we obtain a conformal mapping from the associated 
Riemann surface to the strip ]— K, K[ x IR, which is periodic in the u plane with the 
period 2K' = 2K(1 - m) (see Fig. 

Consequently, we need to look for 2i^'-periodic solutions in the strip ]— K, K[ x M 
with conducting walls. The empty space solution is trivial, 



V{u) = Cu . (42) 

The solution with a point-like charge Q at u = uq (or more precisely at u = UQ+2iK'n 
as required by periodity) is also well known and reads 



V(u) = -Qlog 



0i hr 



AK I 2K 



7T- 



u+Uq— 2K i 



IK 



'2K 



(43) 



where Q\ is a Jacobi elliptic theta function. The result equals the field generated by 
an infinite lattice of positive and negative mirror charges, which are located at the 
zeros of the numerator and the denominator in (43), respectively. Notice that the 
half-period ratio of the solution r = iK' /2K = iK(l — m)/2K(m) is exactly one half 
of the half-period ratio of the mapping p above. 

The above solutions can be used to construct the potential in all cases of interest 
to us. However, the analytic expressions of most physical quantities are quite cum- 
bersome. In some important special cases there are however drastic simplifications. It 
is useful to split the potential into symmetric and antisymmetric parts (with respect 
to the imaginary axis), V = Vs + Va, where 



VM 



V A {u) 



Q 
2 

Q 



log 



log 



01 7T 



u-uq I ■ k' 
iK \ L 2K 



4K 



'IK 



6x ,u_+^-2K HK 



0i hr 



AK 

u—uq 
AK 



2 A 



U—UQ—2K | ■ K 



■IK 



'2KJ_ 



'2K 



u—Uq—2K i • K 
h I 7T % \% — 



AK 



2K 



. u+u*-2K.. K 
01 7T TT? \l 



AK 



2K 



h 



U + U Q 
AK 



'2K 



(44) 



(45) 



Now the symmetric solution may be written (up to a constant factor) in terms of the 
elliptic trigonometric functions, and the asymmetric term simplifies also to 



V s (u) 
V A (u) 



Q 



log 



sc 



U-Uq 



m sc 



u + u* n 



m I nd 



U - Uq 



m) nd 



u + u* Q 



m 



q, »■(- 

— loe 7— 

2 n I U+U 



2K I K 



0i u 



2K 



■Ki 

' K 



(46) 



5 We can also remove the ambiguity of periodicity by mapping the strip to the annulus with the 
exponential map. We, however, find it more convenient to work on the strip. 
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where sc and nd are Jacobi elliptic trigonometric functions and now all the functions 
have the same half period ratio as the conformal map p. The asymmetric potential 
vanishes on the imaginary axis as well as on the conductors. By using the trigonom- 
etry of the elliptic functions, the symmetric term on the z plane may be expressed in 
terms of elementary functions 

V s (z) = V s (u = p{z))\ uq=p{zq) 



Q , 

-2 log 



az\ 1 — 



1 - a 2 z 2 Q - azQ\/l - — Vl - a 2 z 2 



,*2 



x j az\l I - %J\ -a 2 z* 2 + az* Q \jl - —^\-a 2 z 2 



y *2 



v v cr V a" 1 a a 



Vl - a 2 A A -aV 2 



(47) 



The (conjugate of the) electric field and the charge density on the conductors are 
given by 

y s {z) = 



Ps{z) = 



Q\z\ 



a 2 z 2 



[ns (sn-^-siT 1 ^)) + ns (sn-^^+sn" 1 ^))] 



7i\/a 2 z 2 — ly/a 2 — z- 



:Re 



l - a2 ZQ\J0 2 ~ Z 2 Q 



(48) 



where the parameter of the Jacobi functions is m = a -4 and for the charge density z is 
real with 1/a < \z\ < a. Unfortunately, we were not able to express the antisymmetric 
term similarly in terms of elementary functions. 



3 The boundary one-point function 

Consider the boundary one-point partition function: 



Z 0+1 (P,£,N) = Ijj 



' N 

n 



dtj 
2tt' 



n 

l<i<j<N 



iti 



(49) 
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which equals the partition function of Coulomb gas on the circle with an external 
charge £ at w — 1, see, e.g., |6l. The subscript + 1 denotes one boundary charge 



(\w a \ = 1), and zero bulk charges (|w a | < 1)- This case is known to correspond to 
the circular Jacobi ensemble (see, e.g., [16], Sec. 2.9.6; shift U — > U + n), and the 
partition function can be evaluated exactly using the Selberg integral formula. We 
will discuss that in Section |3.2| Since this case is in analytical control, it is a good 



starting point for our scaling limit analysis. 

Recall that we consider the limit of large N and £ with N/£ fixed, described 
by classical electrodynamics. We assume £ > so that the integral (49) is surely 
convergent. Thus the ^-charge repels the unit charges creating a gap in the charge 
distribution near t = where the charge density vanishes. Let us say that the gap is 
created at t = —t c . . . t c . Everywhere else on the unit circle the distribution is such 
that the electric potential is constant. Note that consistency requires that the charges 
at the edge of the gap feel no force. The charge distribution must therefore vanish 
continuously at t — ±i c ; otherwise it would not be an equilibrium situation. 

where a 
finds 



It is satisfying to find out that our analysis reproduces a result of 25 



matrix model was constructed for the circular Jacobi ensemble (49). Ref. 25 



that the spectral measure converges to a probability measure supported by an arc 
of the unit circle, in a limit which corresponds to N, £ — > oo with iV/£ fixed. Our 
calculation can thus be viewed as an alternative, physics motivated, derivation of 
that result of 



25 , and as a test of our approach which can be applied to a variety of 



different cases. 



3.1 Large N calculation 



The electrostatic problem which corresponds to the one-point boundary function was 
already solved above in Section |2.2.2 To illustrate the use of the alternative, more 
general method described in Section |2.2.1[ and to keep the discussion uniform with 



Section |4j we will anyhow outline how the solution is obtained by using the general 
method. 



We will first map the disk to upper half plane, using the function q(w) of (12) 
with C = cot(t c /2): 

tc\ 1 — w 

z i-)- i cot I — ) - = q(w) , (50) 



2 J 1 + w 

see Fig. [T] As the conductor is located at ] — oo, — 1] U [1, oof, we shall choose the 
functions g and h of (13) as 



g(z) = Vl - zVz + 1 ; h{z) = -\fz - W z + 1 



(51) 



It is then straightforward to find the "tilded" solutions described in Section |2.2.1 
We set the charges as 



-N-Z 



at 
at 



z 
z 



(w = 
-i( (w 



v. 



oo) 



(52) 
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where the auxiliary charge £2 is the image under q(w) of the charge at infinity on the 
w plane. In the constraints (23) only the former equation (with the index k — 0) 
gives a non-trivial condition, 

e = 4±L, (53) 

or equivalently sin(t c /2) = £/(-/V+£), whereas the latter one is satisfied automatically 
because of symmetry. By using this constraint together with Eqs. (22) for the above 

arge den 

W^TI (n + o«C 



configuration, the electric field and the charge density on the half plane simplify to 



E z) 



£cv 



Z" 



ir\z\ (z 2 + ( 2 ) 



(54) 



which can be mapped to the disk by the map z = q(w) as shown in (26). 

One can check that the (rather cumbersome) potential matches with the expres- 
sions of Section 2.2.2 both on the half plane and on the disk. However, the potential 
on the conductor U of Eq. ( 29 ) and the "interaction energy" £ of Eq. ( 30 ) reduce 
to very simple forms when ( is eliminated by using ( 53 ) . We find 



where 



U = -[-f(N)-f(N + 20 + 2f(N + 0] 
Zc-t = £[-/(£) + /(2£) + /(iV + £)-/(iV + 2£)] 

f(x) = xlogx . 



(55) 
(56) 



Thus, the total energy reads 
£(ipt) : 



[-2F(0 + F(20 - F(N) - F(N + 20 + 2F(N + £)] , (57) 



where 



F(x) = -xf(x) = -x 2 \ogx 



(5f 



2 " y ' 2 

The resulting leading asymptotics of the partition function for iV — > 00 with iV/£ 
fixed is 

logZ 0+1 (/3,£,A0 



P 



-0£ 

[2F(0 - F(2£) + F(N) + F(N + 20 - 2F(N + £)] (59) 



The charge density becomes, after mapping to the circle, 
P(t) = 




cot" 



cot" 



(60) 
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This result agrees with (5.6) in 25 . Their d is equal to our ratio £/N, which is real 
valued, so that their ^ = 0. The gap angle 0d in |25|, 



sm 



d 



i + d i + {i/N) 



(61) 



is seen to be equal to our gap angle t c . Finally, our charge density (60) can easily 
be rewritten in the form (5.6) of 25 . The result may also be used to calculate the 
next-to-leading correction to log Z. We find 



dtp(t) \og2np(t) 



£log2 + -f(N)+f(N + 



:f(N + 2£) 



(62) 



Therefore, using Eq. (122) from Appendix [Aj the next-to-leading result reads 
logZ 0+1 (/?,£,AO = fL[2F{£)-F(2Z)+F(N) + F{N + 2Z)-2F{N + Z)] 



+ 



2 

2-1 



£ log 2 + ±f(N) + f(N + £)-lf(N + 20 



1 



N + O (log N) .(63) 



Let us make one more comment which is special to /3 = 2. At this temperature, 
retaining the leading term of log Zo+i (the next-to-leading term vanishes), and further 
considering the limit N/ ^ > 1 we find 



J 0+l 



(64) 



This agrees with the leading term of log Zq + \ in the limit of large iV with £ fixed [20] . 
Thus the electrostatic result correctly interpolates between the leading terms in these 
two limits. We have checked this explicitly only for the one-point amplitude, but the 
result is true in general. This is evident since the large iV result of Eq. (64) has 
a simple electrostatic interpretation: the corresponding energy — (£ 2 logiV)/2 is the 
missing self-energy of the small gap created by the external charge in the otherwise 
flat charge distribution [6]. 



3.2 Comparison to the result for arbitrary iV 



The exact result for the partition function (|49j) for general iV and (3 can be found by 
using the Selberg integral [l6l[26] and reads 



7 (R c ,n_TT r ((/3/2)j + ^+l)r((/3/2)(j + !) + !) 

o+i(/u ' } " fi mm + wm + i)] 2 r«/v2) + 1) • (65) 



It is straightforward to check that its expansion at large iV and £ equals Eq. (63). 
For f3 = 2 the result (65) can be expressed in terms of Barnes G functions as 

Zo+l(2 '^' iV) " G(2£ + l) G(N + £ + iy • (66) 
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In particular, either by using the general result of Eq. (65 ) of by using the asymptotics 



of the Barnes G function in Eq. (66), one can check the O (N) term in Eq. (63) indeed 
vanishes at (5 = 2. 



The "exact" counterpart of the continuum charge density (60) (or (38)) is the 



distribution function of unit charges in the Coulomb gas. It can be acquired by 



leaving one of the tj's (say t\) unintegrated in (49). That is 



Pn{U 



Z{£-N)N\ 



' N 

n 

i=2 



dU 
2^ 



N 



i=i 



— e 



n v 

X<i<j<N 



e itj \ 2 



(67) 



where we again fixed (5 = 2. In the large N limit, with £/N fixed, we expect that 

1 



N— >00 _/V 



(68) 



where p{t) was given in Eq. (60). In fact, pjv(^i) (as well as the higher point distri- 



bution functions) is known exactly for (5 = 2 and arbitrary N (see Chapter 4 of |16]). 
By using this result, it is possible to check the limit (68) explicitly. Thus the ex- 



act distribution functions provide yet another way to derive Eq. (60). We have also 



checked the limit by comparing the first few Fourier coefficients of the functions p 



and pn, which are calculable for all £ and N 27 



3.3 Lattice path interpretation 



We first rewrite the partition function (49) as a CUE expectation value, 

Z 0+l (P,tN) = (wi + t*)W(l + e ^)W»\ , (69) 

\»=1 / CUE/3 (N) 

where we have shifted t j — > ti + it. If we set (5 = 2,£ = n, we recognize it as the 



n\ = ri2 = n case of the expectation value appearing in 16 Exercises 8.1 



Z 0+1 (2,n,N) = G l i rd {{lf = 2(.V - ./!}, , v: {Jj = 2{N - j) + r}, v ! 

fe=l,...,2n 

(70) 

with r = 0. The resulting expression Gn^ rd denotes the total number of iV non- 
intersecting single move Id/rd (left diagonal/right diagonal) paths starting from loca- 
tions and ending at lj (see Fig. 3) in 2n steps (segments), in the lock step model 



of vicious walkers 14, 151, a variant of the vicious walker model introduced to model 



wetting and melting. 



4 Higher-point boundary functions 

The boundary two-point partition function can in principle be treated in the same 
manner as the one-point function above. Now we have two charges, £i and £ 2 , on 
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Figure 3: An example of N = 9 non-intersecting left diagonal / right diagonal lattice 
paths of 2n = 10 steps/segments, with endpoints lj shifted from the starting points 
/(o) by r = in the x-direction. 



the unit circle. When these charges are well distinct they form two separate gaps in 
the charge distribution of the unit charges at large N. However, when the charges 
approach each other the gaps will merge at some critical separating angle. Hence, we 
have phases with either one or two gaps in the charge distribution. 

The general solution for the two-point boundary function is complicated. Hence 
we will only discuss two special cases below and restrict to the leading asymptotics. 
We shall also discuss some simple generalizations to higher-point functions. 



4.1 "Antipodal" two-point function 



When the two charges are exactly at antipodal points of the circle, the solution has 
two gaps regardless of the values of £i and £2- We can fix the charge £1 at t\ = 
{yj\ = e ,Tl = 1) and the charge £2 at r 2 = 7r (w 2 = e lT2 = —1). Then the charge 
£1 (£ 2 ) is mapped to Z\ = (z 2 = 00) in the mapping z = q(w) of Eq. (12). For 
a proper choice of the parameter ( (to be solved below) the configuration is that of 
Section 2.2.2 In particular, the system is reflection symmetric with respect to the 



imaginary axis. Hence, the potential on the half plane can be expressed in terms of 



the symmetric solution of Eq. (47) as 



V(z) 



Vs(z) 



=6, 2Q=0 



+ V s (z) 



=£2, ZQ=ioc 



+ Vs(z) 



=-(AT+£i+6), z Q =-iC 



+ V , (71) 



where Vq is a constant. The charge density is constructed by using Eq. (48) and 
becomes 



P{z) 



ix\/a 2 z 2 — l\fc 



(iV + £ 1 + 6)v / i + « 2 C 2 \/^+C^ Sia c 



(72) 
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Requiring it to vanish at z = ±a and at z = ±l/a fixes 



(N + ei)(iV + 6) + 66 + y/N(N + 2^)(N + 2&)(N + 26 + 26) 



If we further fix the constant Vq by requiring that 

V(z) = &logz + 0(l/z) 
as z — > oo, the solution equals the defined in Section 



(73) 



(74) 

It may therefore 

be mapped to the disk by using Eqs. (27). After a straightforward calculation, the 



2.2.1 



potential and the interaction energies (see Eqs. (29) and (30)) become 
U 



f(N + 6 + 6) - 4 [f(N) + f(N + 26) + f(N + 26) + f(N + 26 + 26)] 



(6 + 6 

6 



f(N + 6 + 6) -3/(^ + ^1 + 6) 



|/(iV + 26) 



2 /(iV + 26) - ^F(6 + 6) + ^ [^(26 + 26) + ^(26) + ^(26)] 
-66 log 2, 



(75) 



where again f(x) = xlogx and F(x) = (x 2 loga:)/2. Thus, the leading asymptotics 
of the partition function for the antipodal configuration is 



logZi 



0+2 



(2pt,a) 



2 2 



/3 -F(iV + 6 + 6) + i ? (6 + 6) 



+j [F(iV + 26 + 26) + F(N + 26) + F(N + 26) + F(N)] 



4 



[F(26 + 26) + F(26) + F(26)] + 66 log 2 ^ . 



(76) 



4.2 Two-point function with equal charges 

Let us then discuss a slightly more complicated special case where the two charges 
are equal 6 = 6 = £ Di rt their locations are arbitrary. When the angle r = r 2 — T\ 
between the charges is smaller than some critical angle r c , there is only one gap. 
By rotational symmetry the two charges can be taken to lie at e ±JT//2 on the unit 



circle. Then the mapping (12) with ( = cotr/4 maps the 6 criar g es to z = ±1 
and the w = 00 point to z — —i( = — zcot(r/4). We can then use the methods of 
Section EXH with 

g(z) = Vz - b\/z + b , (77) 
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where 1 > b > and we used the fact that the solution must be symmetric with 
respect to the imaginary axis. The constraint equation for the parameter b simplifies 
to 

2£ iV + 2£ 



(7f 



Following the procedure described above and in Section 2.2.1 and after a lengthy 
calculation, the energy can be expressed as 



£(lg) 



±F(N)+F(N + 2£)-±F(N + 4£) 



+ ^N{N + 20 log (C 2 + 1) + 6£ 2 log 2 - log (2 sin 



(79) 



For t > t c there are two gaps. Then we can use the same choice for the mapping 
z = q(w) to the half plane as in the one gap phase, with ( = cot J. Now we choose 
in Section \2. 2. II 

g(z) = yfz + ay/ z — ay/ z — b\J —z — b , (80) 
where we used symmetry to fix two of the parameters defining the end point locations 



and a > 1 > b > 0. The remaining parameters a, b are fixed from the equations (18) 
which simplify to 

(N + 20 C V^iVT^¥ = 2£y/eT^y/eT¥ 

C 2 (a 2 - 1) II (1 - b 2 1 1 - b 2 /a 2 ) + (a 2 + C 2 ) Ren (l + 6 2 /C 2 1 1 - b 2 /a 2 ) =0 (81) 

where IT(n|m) is the complete elliptic integral of the third kind. Taking the real part 
of II in the latter equation removes the branch cut ambiguity, which appears since 
the argument n — 1 + b 2 /( 2 of U(n\m) is larger than one. 



The critical angle r c may be extracted from (81 ). For r = r c either of the conduc- 



tors should disappear and the gaps should merge together, i.e., we need to study the 
limit b — >■ or a — )■ 00, which yields 




(iV + 20 2 tan 2 (f) +40 



N(N + 4£) 



N 2 - N(N + 40 cos (f ) + 4^ + 80 



(82) 



After a tedious calculation by using the formulas in Section 2.2.1, the energy can 
be expressed as 



£(2g) 



2(c 2 + ik 2 t , (c 2 + i)^+2o 



y/l^S 2 



-Jl + 



(83) 
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where 



lim 

e4.o 



a , Va 2 - z 2 Vz 2 - b 2 V^TVT^¥, 
dz 7Z2 TY7Z2~i 7* + on , ^ lo § e 



Z 2_1)( Z 2 + C 2) 2 (1 + C 2 ) 

y/a* - Wi - b 2 log ( ^fey^iy ) (6 2 + C 2 ) ^(sin- 1 ^ 



a 2 -l 



2_ fe 2 



II - 



S) 



2(C 2 + 1) 
+ C 2 )n(^;sin- 



«(c 2 + i; 



a 2 -l 
a 2 -fc 2 



II - 



6 2 



'i-& 2 m i 



a(C 2 + l) 



-4; sin 1 



a 1 -! 
a?-b 2 



b 2 



lim 

40 



a(C 2 + l) 

^ I 9 , -iW 9 A9\ T n/-//-9 , i\ ^°S e 



( y 2 + 1) (y 3 _ C 2) 2C(C 2 + 1) 

v^Tev^Te log ( jg^ y ) , (6 2 + C 2 )F(cot-(aii-S) 
2C(C 2 + i) 

(l-& 2 )n(l-i; cot- -| 



+ 



a(C 2 + l) 



a (C 2 + 1) 

(& 2 + C 2 )n(^g^;cot- 1 (i)|i-g 



log 



2£ 2 log 

— £ 2 log ( 2 sin 



2\/2C 
C 2 + l 
r 
2 



a(C 2 + l) 

-^(iV + 20 2 log(2C) 



The calculation was done by using integral tables and checked by Mathematica. 
F(u\n) [II(n; u\m)] is the incomplete elliptic integral of the first [third] kindj^] 
The one and two-gap solutions are combined to give the total energy as 



£(2pt) = £(lg)0 _ r) + ^(2g)^ (r _ r j 



J5) 



in Fig. |4j Notice that the solutions join smoothly at r = r c . When the interaction 
term £g is included, the r dependence is relatively flat for large separations (r close 

to 7r). 

Let us then compare the leading term of the partition function (at £i = £2 and 

= 2), ' 

Z 0+2 (£,iV)~exp(-2£( 2pt )) (86) 

to the exact result and earlier approximations. The exact partition may be calculated 
at small N by noticing that the partition function at f3 = 2 defines the Toeplitz 



6 The result could possibly be simplified a lot, if the constraint equations could be solved and 
substituted into (84). Such simplifications take place for simpler configurations. 
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Figure 4: The boundary two-point total energy as a function of r. We used £ = 6 and 
iV = 16. In the left hand plot the energy of the mutual interaction of the pointlike 
charges 8% was excluded. The critical angle is given by r c ~ 1.157. 



determinant (see [9]' 



Z 0+2 (f i, 6; AO = |e* T1 - e^l 266 det 7>(/i) (87) 



for the function 

hit) = | e lt - e iTl | 26 | e 4 * - e iT2 | 26 . (88) 

The partition function can be evaluated in the limit iV — > oo with £ fixed by using 
the asymptotics of Toeplitz determinants 20 . Leading term given in the literature 
can be written in terms of the Barnes G function as 

z (t r - n) ~ IT G ^ a + 1)2 G{N + 2ia + l)G{N + 1} (89) 



when applied to our case. Based on physical considerations we suggested [6] the 
improved formula 

, r(iv + ?1 + 1 )r(iv + ?2 + 1 ) 



T(N + 1)T(N + Ci + 6 + l)G(iV + 6 + 6 + I) 4 

x n rlt + +?\ G ( N + & + E & + + e. + 1) ■ (90) 

fJi G{2^ a + 1) V 



Notice that both above approximations are independent of r. They are expected to 
work well when the charges are close to antipodal even for relatively large £, but 
they will break down for r < £/N (see the discussion in (6j), as we shall see. In the 
region of small r we may in fact approximate the partition function by separating the 
interactions of the external charges and then joining the charges into a single one: 

Zo+2(Zi,&N) ~ le^-e^f^Zo+^ + ^iV) (91) 



,266 ^(Ea & + I)' ^(iV + 2 E Q ^ + l)G(iV + 1) 

1 G(2L( a + l) G(iV + E a ^ + l) 2 
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Figure 5: Comparison of the electrostatic large £, N approximation of Z0+2 to asymp- 
totics of the Toeplitz determinant and to the exact partition function at /3 = 2 as a 
function of r. We used £ = 6 and N = 16 (N = 16 and iV = £) on the left (right) 
hand plot. The solid black curve is the leading term obtained by electrostatics given 
in Eq. (86). The exact result is shown by the dashed blue curve which is hardly 



visible since it overlaps with the leading term. The gray dash-dotted horizontal line 



arises from the asymptotics of the Toeplitz determinant, Eq. (89). The horizontal red 



dotted and the green dotted lines are the improved asymptotics of Eq. (90) and the 



small r approximation of Eq. (91). 



where we used the exact result for the one-point function on the second line (or 
equivalently either of the approximation schemes above, since they all give the same 
result for the one-point function). 

All these approximations are compared in Fig. [5j As expected, the asymptotic 
formulas obtained from the asymptotics of Toeplitz determinants work well for the 
region of large r where the partition function is almost constant. Actually, the im- 



proved asymptotic formula (90) (dotted red line) gives a much better approximation 



than the standard leading term of the Toeplitz asymptotics (89) (dash-dotted gray 



line). The leading term at large N (86) (black curve) interpolates between the other 
approximations, and follows closely the exact result (blue dashed curve), which is 



obtained numerically from Eq. (87). 



4.3 Symmetric n-point functions 

Let us consider a symmetric configuration (see Fig. [6]) where n external particles have 
all equal charges £ and are located at angles r a = 2ira/n, where a = 1, . . . ,n. For 
a continuum charge equally distributed into n symmetric blobs between the external 
charges, the potential can be immediately written in terms of the boundary one-point 
solution above. Namely, 

V inpt \w) = V {lpt \w n ) 
pto*)(t) = np (lpt) {nt) . (92) 

We show the potential and the electric field for a choice of parameters in Fig. |6j 
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Figure 6: The potential Rey(" pt '(w) (left) and the absolute value of the electric field 
(right) for n = 7 with N/C, = 6. The continuous charge distribution is denoted by the 
black arcs. 



The value U of the potential at the conducting planes is unchanged, as well as the 
potential ReV^(u; a ) felt by the external charges, which fixes £ c -£- Notice that while 
conformal mappings conserve charge locally, in«ji-> w n the total charge is multiplied 
by n since there the mapping takes the charges to n distinct images. Hence, we 
actually should start by a one-point solution with jV( lpt ) = N/n = N. Consequently, 
the energy can be written as 



g(npt) 



n 



2F(f) + F(2£) - F(N) - F(N + 2f) + 2F(N + £) + ^ , (93) 



where the term arising from the mutual interactions of the external charges reads 



n-l 



^log|e 2 ™ /n -l| 



< 2 



logn . 



(94) 



o=l 



Up to this additional term the result agrees with the approximation suggested in [6] , 
which was based on the limit iV — > oo with ^ fixed. 

It is clear that the mapping w t— > w n can also be applied to more complicated 
solutions than the boundary one-point function. For example, we can map the an- 
tipodal configuration discussed above by, say, w H- w 2 , to obtain the asymptotics for 
a (quasi)-symmetric four-point partition function with two equal charges £i at r = 
and at r = it, and charges £2 at r = n/2 and at r = 3n/2. The total energy for this 
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system becomes 

s = 2s^\ N ^ N/2 -(e l+ e 2 )io g 2 

= 2F(JV/2 + £i+6)-2F(£i+6) 

~ [F(N/2 + 2Ci + 26) + F(N/2 + 2£i) + F(JV/2 + 2£ 2 ) + F(iV/2)] 

+ ^ [F(26 + 2^) + F(2&) + F(2&)] - (6 + 6) 2 log2 , (95) 

where the interactions of the external charges within themselves were treated sepa- 
rately. 



5 The bulk one-point function 

We shall now point out a connection between the results of our previous article [7] 
and the Poisson kernel. Consider the bulk one-point partition function: 



Z 1+0 {w,£;N) 




' N A+ 

naii 



i=i 



2vr 



n y 

l<i<j<N 



we 



-iti\Pt 



(96) 



CUEp(N) 



which equals the partition function of the Coulomb gas on the circle with an external 
charge £ at a fixed complex position w^j When £ = — (N — l) — 2//3, the integrand is 
a special case of the Poisson kernel, 



1_ / det(l-gtg) 
( } C Vl(det(l - S-tS 1 ))! 2 



(97) 



when the matrix average is proportional to the unit matrix, S = wIn, and when 
the random matrices S belong to the CUE. In condensed matter physics, this case 
is used as a model for transport across a disordered cavity, connected to the outside 



by non-ideal leads [28] (see 29 for a general review). Here N denotes the number of 
scattering channels. In this special case, the probability distribution function becomes 



N 



P({ti}) 



dt 



Z 1+0 N\ Al 2 



(9S 



l<i<j<N 



where we recognize the normalization factor Zi + q as (96). In this case, when (3 = 2, 
£ = —N is a negative integer, < w = r < 1, (96) has been calculate^ in 30 , with 
the simple result 



'1+0 



>,£;AT) 



l-r 1 



(99) 



7 Notice that by rotational symmetry Zi + $ is independent of argw. 

8 Note that |30| assumed N > |£| but this condition is automatically satisfied with the above 
choice of £. 
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In the case of positive £, (96 ) is much more complicated to calculate. In the double 



scaling limit we have obtained the result [7|, 

Z 1+0 (r,£;N) 



-2£+V(N) 



where the exponent is 



£ 



X 2 



Q( r c - r)y log (l 



r 2 ) 



(100) 



(101) 



+9{r — r c ) 



{N + 2if 1 + x 
-log 



4 



1 + 6(r) 



N\ l-x 

X log i3^y 



4x 



H log 

2 S (l + 5( r )) 2 J 



where 



N 



6(r) 



1 — r 



(102) 



iV + 2£ ' w 1 + r ' ^ iV + £ ' 

and the term T>(N) includes contributions from subleading terms, in particular from 
terms of O (N°). The function in the exponent has two possible forms, depending on 
the location of the bulk external charge. If it is sufficiently close to the unit circle 
(r — > 1), a gap [— t c , t c ] opens up in the continuous charge density p(t), with the size 
2t c given by 



tr 



COS 



5(r) 



(103) 



where x — + Thus t c vanishes at r = r c which signals the transition between 



gapped and ungapped charge distributions, reflected in the form of (101). 



The result (100) for positive £ can be applied to calculate the cumulative distri- 



bution function for long lattice paths in the discrete polynuclear growth model, as we 
will discuss next. 



6 Bulk n-point functions and polynuclear growth 

Next we study the gas with n external charges in the bulk. Recall that the partition 
function ( [49] ) becomes the CUE expectation value (we repeat the equation ^ for 
convenience), 

I n N \ 

z n (J3; {U; N) = J] K - ™ b \ mUb ( II II I 1 - wae-^l*" ) • (10 4 ) 

a<b \a=l i=l / CUE^j(V) 



We point out that (104) is related to the weight of a lattice path in the discrete 



polynuclear growth model. We will follow the presentation of the model in 17 or the 
review [16] where more details can be found. Here we just summarize the basic idea. 
One starts from an ensemble ofmxm weighted non-negative integer matrices X = 
(xij). Each entry of the matrix carries a weight factor {1 — aibj)(aibj) Xli , interpreted as 
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coming from a geometric distribution (hence the normalization factor (1 — cubj)). One 
then rotates the matrix counterclockwise by 45°, and starts reading the entries from 
the bottom corner, and growing a box diagram following a set of rules. Basically one 
starts from the y = ground level and keeps adding rows of vertical boxes, the height 
of each box given by an entry of the matrix. The addition of a box is a "nucleation 
event". The matrix entries are read one horizontal row at a time starting from the 
bottom corner. At each step, the previously constructed part of the diagram grows 
wider by one step to both left and right, before the new boxes are added on the top. 
Thus the box diagram will grow into a pyramid shape. The left and right sides of 
the boxes carry weights b Xij and a Xij , respectively. As the boxes grow wider, they will 
sometimes overlap, in which case the overlapping part of the diagram is moved to a 
new sub-diagram below the ground level, y < 0. In the end, when the box diagram is 
complete, the left and right edges of the profiles of the main diagram and subdiagrams 
are interpreted as pairs of weighted non-intersecting lattice paths of 2m — 1 steps. 
The profiles with endpoints at y — — I + 1, I — 1, . . . , m are called level-/ paths. See 
Fig. [7] for a sketch of the process (following [l6]). 
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Figure 7: A sketch of mapping a non-negative integer matrix to a box diagram ac- 
cording to the rules of the discrete polynuclear growth (compare with [i~6]). 

An interesting random variable in the polynuclear growth model is the maximum 
height of the level- 1 path (from the profile of the main part of the diagram, starting 
from and ending at the ground level y = 0), denoted h u . A fundamental statisti- 
cal quantity is its cumulative probability density Pr(/i n < N). This is where the 
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connection to circular unitary ensemble can be established. One can show that 

m I m N \ 

<n)= - ai&i) • ( n ri( i - - b ^ ) ■ ( ios ) 

i,3 \j=lk=l I cue^at) 

Now consider the special case of left-right symmetric weighting, bj = aj. Furthermore, 
suppose that m = rti\ + m 2 + • • • + m n for some n, and make the following special 
choice for the weight factors, 



imi — ' I; Q"mi+1 



'mi+m2 — ^2) • • • j Q"m— m n +l 



Then (105) becomes 



n AT 



Pr(/> n <iv) = n( 1 - r » r >r° 

a,6 



m 6 



ruin 

\a=l k=l 



-it k \2m a 



r n ■ 
(106) 



(107) 



CUE 2 (Af) 



We recognize this as a special case of (104) with 



(3 = 2 ; £ a = m a ; w a 



1. 



n 



(108 



adjusting the normalization factors. Thus, the bulk n-point functions compute the 
cumulative probability densities of the level- 1 paths in the discrete PNG model. Note 
that a bulk two-point function with r\ = 0, r 2 = r formally becomes equal to the bulk 
one-point function, which we discussed in the previous section. In this case PNG 
becomes independent of £i = mi, and denoting £ 2 = m 2 = m we find 



A' 



Pr(h n < N) = (1 



r 



2\m? 



\k=l 



— re 



(l-r 2 ) m Z 1+0 (r,m;N) 



CUE 2 (A0 



(109) 



The expectation value is the one-point partition function (96) with w = r, /3 = 2,£ 



m. In the double scaling limit N, m 3> 1 with N/m fixed, we derived and approximate 



result (100), so the cumulative distribution function becomes 

Pr(h n < N) = (1 



r 2 ) m2 e- 2£ 



;no) 



with the exponent (101 ) and up to the subleading corrections V(N) in (100) which we 
shall neglect. We expect that form to be valid even if we keep m fixed and vary N, as 
long as N, m ^> 1. In the growth model, the weight factor r has a fixed chosen value. 
On the other hand, for the cumulative distribution function, we need to establish that 
Pr(h n < N) — > 1 as iV — > oo. Since r is fixed, at some point it will become less than 
the critical value r c = N / (N + 2m), when N increases. When r < r c , the exponent S 
becomes 



£ = ^log(l-r 2 ) 



111) 



29 



which is independent of N. We interpret this as a sign that the cumulative probability 



distribution has saturated, indeed (110) becomes equal to one (up to higher order 
corrections). Thus Pr(/i n < N) — > 1 as N — > N c , where 



Nr. 



2mr 



112) 



found from (see (102)) 



(113) 

N c + 2m K J 

This means that N c is a probabilistic upper bound for the height of level- 1 paths, in 
the sense that cumulative distribution function saturates. For N < N c , the cumulative 
distribution function is 



Pr(/i n < N) 



X 



i + x 

1 + S(r 

N 2 /2 



1 - 5(r) 



(N+2m) 2 /2 



4x 



;ii4) 



Fig. |8] shows an example plot for m = 10, r — 0.5 (so that N c = 20) in linear and 
logarithmic scales. We also show the exact probabilities which are easy to obtain from 
Eq. ( |109[ ) for not too large iV by evaluating numerically the corresponding Toeplitz 
determinant 

Z 1+0 (r,m;N) = detT N (h) 
-it\2m [ com p are t Eq. (87) . 



;ii5) 



for the function h(t) 



re 




Figure 8: The cumulative probability distribution function (114) in the scaling limit 
(solid line) compared to the exact result (crosses). 



6.1 Lattice path interpretation 

The bulk n-point functions at = 2 and with charges on the real line also admit a 
lattice path interpretation which is very similar to that of the boundary one-point 
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function discussed in Section [3~3 Actually, we may again write a connection to iV non- 
intersecting single move Id/rd (left diagonal/right diagonal) paths with 2m = 2 J2 a m « 
steps with the same starting and ending locations as for the boundary function (see 
Fig.pl: 



n 



.a<b 
Qld/rd 



r a - r b \ 

(o) 



-2m a m b 



N n 



'2m 



Z„+ ({r a }, (».}: N) = { Yl Yl |1 - r,t 

k=l a=l 

2(N - ./)}, , y: {lj = 2(N - ./)}• , v) , 



-it k \2m a 



CUE 2 (iV) 

(116) 



but now we need to add some nontrivial weights to the paths (see 16 , proposition 
8.13). We denote the weights of the left (right) diagonal paths at step k by w^. The 
weights for the m first steps are 



w-. 



w 



mi+l 



W mi = r l 



10 



mi+m.2 



r 2 



w 



m—m n +l 



u ■■»■. 



while for the remaining m steps we flip + <-)•—, i.e. 



w 



T 

k+m 



W 



k i 



, m 



(117) 



(118) 



According to (116) the n-point partition function is then proportional to the total 



weight of all N non-intersecting Id/rd paths with above boundary conditions and 
weights. Letting all r n — > 1 we recover the result for the boundary one-point function 
(with charge m) as all the weights become equal to one. Notice also that for the n = 1 
case we can immediately extract the leading behavior of the total weight G^m^ m the 
double scaling limit by using Eqs. (100) and (101). It would be interesting to know 



what other interpretations of the bulk n-point function (96) exist for other values of 
(3^2 and noninteger £ a . 
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A Partition function to next-to-leading order 



Consider the two-dimensional Coulomb gas on the circle, subject to an external po- 
tential £V ext (z). The system is defined by the Hamiltonian 

N N 

h = -J2 lo § I e * - eltj I + £ E v ^ ( elh ) ' 19 ) 

i,j = l i=l 

where the positions of the unit charges are Zi = e lti . The partition function is then 

■ j i 

We consider the limit where N, £ — > oo with N/ £ and V ex t fixed. This is a gener- 
alization of the limit which was discussed in the main text to an arbitrary external 
potential and general (3. The aim of this Appendix is to sketch how to calculate the 
first two coefficients in the series ^ where the coefficients depend on £/N and may 
also include logarithms of N. In particular, we show how to relate Co to the electro- 
static energy in the continuum limit. We shall do a direct brute- force calculation of 
the coefficient C±, which is shown to have a simple structure. A more sophisticated 
method could use an effective theory of continuous charge distributions |31~] . 

Recalling the saddle-point approach (J5|, we write first log Z(/3, £, N) = —(3H + 
Af, where Hq is the (global) minimum value of the Hamiltoniarj^] and A/ is the 
contribution from fluctuations around the minima. For nonzero external potential, it 
is hard to calculate H explicitly. Hence we introduce a further correction term by 
writing —{3H Q = —(38 + Ad, where £ is the electrostatic energy in the continuum limit, 
and Arf is a discretization correction. The continuum energy — (3£ scales as N 2 and 
gives exactly the leading term CqN 2 , while the next-to-leading coefficient C\N can be 
extracted from the various correction terms A. We shall calculate the discretization 
correction at O (N) by further dividing it to two separate contributions: one arising 
from "self-energies" of the point charges, and one from differences in near-neighbor 
interaction strengths. Hence we have, at O (N), 

log Z(P,£,N) = -P£ + A f + A sc + A nn . (121) 

For more precise definitions see the analysis below, where we discuss the various terms 
one by one. 

Notice also the equilibrium charge densities of the continuous and discrete systems 
may be different at next-to- leading order in 1/N. One could naively expect such a 
deformation to cause an additional O (N) correction term. However, recall that £ is 
minimized by the leading order charge density of the continuous system. Therefore 

9 This step may be subtle if the external potential is complicated, and H has several separate 
minimum configurations. This can occur if the support of the charge density is disconnected in the 
large N limit. One can check that our analysis remains valid when global minimum configurations 
are used both in the discrete and continuous systems. Actually, we will not use the results of this 
Appendix in cases where this subtlety appears. 
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the next-to-leading order change of the total energy due to any perturbation of the 
charge density vanishes. 

In summary, using the results for the correction term A from Eqs. (137), (141), 
and (144) below, we obtain the leading term CqN 2 = —08 as in ([6]), followed by 



C X N = A / + A sc + A nn + O(iV ) 

6—2 f 

' dtp(t) log27r p(t) + 



2-/3/3, 8 

— + 2 l0g 2- l0gr 







+ 1 



N ,(1221 



where the total energy of the continuum system in equilibrium S and the corre- 
sponding charge density p(t) can in principle be calculated by using the methods of 
Section [2.2.1 Notice that C\ vanishes for = 2. 



A.l Dyson gas 

To motivate the calculation of Ay below, we consider the case V eyit = 0, i.e., the Dyson 
gas |8j. Even though we will be interested in the limit of large N with fixed, it is 
illustrative to consider the low temperature behavior (0 — > oo) first. The partition 
function is given by 

w-w = M/n|nk"*-^r p*» 

i Kj 

In this case Cq = (consistently with the expectation from electrostatics). For large 
0, the partition function may be calculated by using a saddle-point approach. By 
rotational symmetry, we fix t\ = 0. There are (N — 1)! equilibrium configurations, 
given by the permutations of the configuration tj = 2ir(j — 1)/N with t\ fixed. The 
value of the integrand at the equilibrium is given by 

^ | e 27riCj-i)/iV _ e 2ni(k-i)/N^ _ _ ^ 124 ^ 
j<k 

Because of symmetry we may fix = t\ < ti < ■ ■ ■ < tjv < 2-n and expand around 
one of the minima. By using 

^ w (125) 

we can write 

1 f -A- dsj 



Z(0,N) = ^ J v H^e-^» (126) 
= Ij\T^ J JJ^e-f^-i^fl + O^JV 3 ^)] , (127) 
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where Sj = tj — 2ir(j — l)/N are the deviations of tj from the equilibrium configuration, 
AH is the Hamiltonian minus its value at equilibrium, the matrix M is understood 
to be evaluated at the equilibrium point, and the region of integration V is fixed by 
< t 2 < ■ ■ ■ < t N < 2tt. 

In the end we will consider the behavior of Z at large N. Therefore, we will 



analyze the dependence of the various terms in Eq. ( 127 ) on N in addition to their f3 
dependence. For large N the dominant contribution to the matrix M comes from the 
charges within distances \e ltj — e ltk \ ~ 1/N. Therefore, the (near diagonal) elements 



of the matrix Mjk behave as Mjk ~ N 2 . The integral (126) is dominated by the 
region where the argument of the exponential function is 0(1). Thus, for large fi 
(and N), the variables scale as Sj ~ l/(y/]3N). The boundaries of V lie at Sj ~ 1/N. 
Therefore, for large the saddle-point approximation works, and we find 



Z(p,N) = N^ 



N(27T0)( N - i y 2 VdetM 



1 + 



'1281 



Notice that the O f^fj correction term vanishes in the St integration since the in- 
tegral is odd. The O (j^j term might have a sizeable (polynomially) iV dependent 
coefficient. 

If we take N — » oo the determinant det M may also be evaluated. In this limit M 
becomes a Toeplitz matrix: 

„ f - 2 C(2); j = k 



W M]k ^\ j^k with (j-fc) 2 «iV 2 . ( 129 ) 



More precisely, the elements near the diagonal, which give the dominant contribution 
to the asymptotics, approach the Toeplitz formj^] The asymptotics of det M may be 
solved by using Szego's limit theorem for Toeplitz determinants. We get 

" AT \ 2N 

_ ) JV(log27r 2 -2)+0(logiV) 



detM = l^-J e ^{^-2 )+ ^ og i V ^ (13Q) 
and consequently, taking first — > oo and then iV — > oo, 

log Z((3, N) = ^N log JV + N Q log - l - log tt + 1 J + O ($° log N)+o(j\ , 

(131) 

where the O (jjjj term depends polynomially on N. Now, let us compare this to the 

10 Notice that since the charges lie on a circle, the distance \ j — k\ needs to be understood modulo 
the matrix dimension, i.e., the elements near the upper right and lower left corners of the matrix 
are also relevant. This is also required for M to be asymptotically of the Toeplitz form. 
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exact result 



log Z((3,N) = log 



r(f + 1) 

N 



N\T (f + 1) 
8 — 2 

N\ogN + 



2 
1 



P 



^(l0g27T-l) + ^l0g^ 



logri £+i 



+ 2 l0g 2 + ° 



N 



(132) 



For large N and 8 one may check that (|13 1|) is indeed correct. The correction terms 
are found to be O (N°) + O 



j J so that the O (logiV) terms actually cancel. 

Let us then discuss the limit N — > oo with fixed (3. In this limit the saddle-point 
approximation breaks down. However, we may extract the result for the integral in 



Eq. (126) in this limit by comparing to (132) 



N 



ds 



, — 2tt 

V j=2 



i e -pAH({ Sj }) 



n 

exp \ -iVlogAM 



N + O (N°) 1(133) 



2 ^ + ^.o g ' 3 -io g r^ + i 

/3 2 & 2 & V 2 

This result will be useful in the calculation for general V ex t below. 

Let us conclude the saddle-point analysis of the partition function by an important 
observation. From above calculation we learn that only variables Sj, sj~ with \j — k\ <^ 
N were coupled in the saddle-point integration, as seen from the structure of M. 
That is, only the near-neighbor interactions are relevant when fluctuations around 
the extremal configurations are considered at large N. While we were not able to 
calculate the integral for general j3, this statement remains true independently of /3, 
and also in the presence of an external potential, as we shall sketch below. 



A. 2 Fluctuation corrections 

The treatment of the term A / is the most tricky one, and we will only sketch how to 



obtain the result. In analogue to (126) we may write 

r N rls 

Z(J3, £, N) = e-? H ° / FT ^l e -^ms 3 }) ; (134) 

where AH is again the difference between the Hamiltonian and its equilibrium value, 
and V is defined by t\ < t 2 < ■ ■ ■ < £jv < ti + 2n. We assume that V^t breaks 
the rotational symmetry, and hence fixing t\ is neither possible nor necessary. For 
sufficiently large N, the equilibrium condition can be described by a continuous charge 



11 Interestingly, by doing the comparison instead in the limit j3 — > oo with N fixed we find the 
exact result detM = 2T(N) 2 / (N2 N ). 
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distribution. Let us denote the density by p and assume normalization f ^p(t) = 1. 
Then p = 1 in the absence of external potential. Let us do the change of variables 



Uj = p{tj)Sj = PjSj 



We obtain 



AT 

n 



^ e ^AH({ Uj }) 



2tt 



(135) 



(136) 



The claim is now, that the remaining integral becomes If of (133) for large N. We do 



not have a precise proof (except for large f3 where the saddle-point method works), 
but one can understand why the result arises. As noted above for V cxt = 0, only near- 
neighbor interactions matter when fluctuations around the minima of the Hamiltonian 
are considered. Locally, within distances 0(1/N), in the minimum configuration 
the distances of adjacent charges are to first approximation constant, and equal to 



2ty /(Np). Hence the substitution (135) removes dependence of the Hamiltonian on 
p for interactions within distances O (1/N). At O (N) in Af we may replace AH in 
(136) by its Dyson gas counterpart in If of Eq. (133), and obtain 



Ai 



^logpj+logl) 



dt p{t) log 2-Kp{t) + 



2-/3/3, P . V (P,« 

log log 1 hi 

/3 2 s 2 s I 2 



N , (137) 



where the term involving the charge density could be written as an integral over the 
continuum charge density up to higher order corrections. Notice that the dependence 
on the external potential appears only through the modification of the charge density. 



The above arguments which lead to the result (137) may be illustrated by consid- 



ering the first term in the saddle-point calculation, which is the leading one for large 
f3. The matrix is found by taking the second derivative of the Hamiltonian and 
evaluating it at the equilibrium. The external potential appears explicitly only in the 
diagonal elements, which read 



l 



Me I 



d 2 



■Xtj 



(138) 



The ^-dependent term is suppressed by 1/N with respect to the near-neighbor con- 
tributions, as expected, and can be neglected. Therefore the approximate result for 
Mjj~, arising from near-neighbor terms, can be written in terms of p as 



4vr2 
N 2 



M. 



-2((2)p 2 + O (i) ; 



PjPk 



j = k 



with (j — k) 2 iV 2 



(139) 



At leading order we were able to use the geometric mean density y/ pjpk for the near 
diagonal terms. By making the substitution (135) in the saddle-point integral 



we 



confirm that dependence on p appears only through the change in the integration 



measure. Hence the saddle-point result is in agreement with (137) 
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A. 3 Self-energy corrections 



The terms A so and A nn arise from the difference of the minimum energies of the 
continuum and discrete systems. We calculate them by taking the continuum system, 
by compressing suitable clumps of the continuous charge into point charges, and by 
studying the change in energy. The easiest term to calculate is A se , which arises from 
the self energies of the clumps that are absent in the discrete system. The self-energy 
of a single clump is given to first approximation by 

1 f s/2 3 1 

-8/2 



Ese= 2JP dxdylog \x - y\ = -- + - logd , (140) 



where 5 = 27r/(Np) is the size of the clump. Summing over the clumps gives the 
O (iV) result 

A se = E sc , = ^ log N - ^ log 2vr + + | £ log p 3 . (141) 

3=1 0=1 

A. 4 Near-neighbor interactions 

The term A nn may be calculated similarly as the self-energy corrections. The change 
in energy between two clumps is 



1 f s/2 

AE nn = logd- — dxdy\og\d + x-y\ (142) 

J -5/2 



fS/2 

-6/2 

where d > 5 is the distance between the clumps. For d ~ 5 ~ l/N, we find 

AE nn = l-h m - l) 2 log f 1 - —\ + 2m log (l + -) (143) 
22 \ m 2 J \ m J 

with m = d/5. Summing over all pairs of clumps gives 

A nn = (3NJ2 ^nn = ^ log 2tt - . (144) 

m=l 

It is not hard to see that contributions for d ^> l/N are suppressed. This is signaled 
by the convergence of the sum over m in the above calculation: the region m > N 
where our approximation for AE nn breaks down does not contribute at leading order. 
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